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1. Introduction 


Further progress in industrial electrochemical generators’ designs utilisimg SOFC considerably depends on 
problem of obtaining fuel cell materials with predefined properties [1,2]. Due to this fact, considerable efforts of 
designers are directed on creation of optimum composition «anode-electrolyte-cathode», and on microconstruction 
of SOFC materials (selection of phase and element composition, defect structure formation, etc.). In the second 
problem nanocomposite systems having unique properties are of greatest interest [3]. Utilising of such materials as 
solid electrolyte may allow significantly improve cell's parameters, in particular, increase ionic conductivity [4,5]. 
However mechanisms leading to this effect were not determined singly. 

In this article the influence of nanoparticles ensemble with work of exit that different from the rest of the matrix 
on the change of electron density in SOFC materials have been examined theoretically. The calculations of 
concentration distribution of free carriers around single particle and for microparticles ensemble were fulfilled. In 


conclusion possible directions of forming of nanostructures conformably to SOFC materials discussed. 


2. The Electron Density Distribution Around Spherical Particle. 


Consider the single spherical particle of radius ro in nonmetallic material. In general case the difference of 
thermodynamic (count from Fermi level) works of exit of the matrix Ymat and of the particle par Qo= Xpar - X% 
mat # O 

If @o > O, energetic zones in matrix near the contact "distort up", and around particle the impoverished by 
carrier's area arise. If 0 <0, zones "distort down", and the enriched by free carrier's area appear. 

Evaluate quantitative parameters that system. Introduce the magnitude @ = - eU (e - electronic charge, U - 
electrostatic potential), which has sense of potential energy of electron in electrostatic field created in matrix by 


particle. Call the magnitude @ by potential that submit to equation [6]: 
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with boundary condition: 
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at the boundary S "matrix - particle". 
In equations (1, 2) no - the density of carriers far from particle, ¢- matrix permittivity. 


The density of carriers around particle determines by known expression [6]: 
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In flat case the analytic solution (1) is known [6] and describe widely used Shottky barrier. 


In considering spherical symmetric case the equation (1) has the form: 
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and boundary condition (2): 
lM), = %0| (5) 


For positive and negative of the solutions of equation (4) will be different. Consider separately these cases. 


a) Matrix impoverishment (@ off > 0). 


As far as all physically interesting cases suppose Ch kT, we can neglect by last term in right part of equation 


(4). 
For determination of second boundary condition in the equation (4) introduce spherical «influence zone» of 


given particle of radius R. Lets accept that electric field at the boundary «influence zone» equals zero: 


= (6) 


For numerical estimations of obtained results we need to determine parameters, incoming in equations (4) - (6). 
Doesn't concretizing material, for calculations we will use middle parameters. 
The solution of equation (4) with boundary conditions (5), (6) have the form: 
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Substituting the finding dependence = @(r) in equation (3), we receive the distribution of electron density in 
space around single particle. For illustration the found dependencies are shown at fig. 1 the calculation n = n(r) 


were made for value @o= 0,3 eV and different values ro. 
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Figure |. Density of free carriers vs. relative distance for different dispersion particles sizes (T = 1500 K, @o= 
0,3 eV, Q=1 %, € = 10, no= 10! cm”). 
r-Ip 
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to see, that ultradisperse particles are of the greatest influence on concentrations' changes (compare Fig. 1). As was 


For comparison the significance n has been given in dependence from relative distance & = It is easy 


shown by calculations «effect of impoverishment» will increase when @o increases under other equal conditions. 
Parameters used for calculations: T = 800,1200,1500 K; no = 10!°...10!® cm; ¢ = 10; ro =10...2500 nm; -@o 
(enrichment) =0,2...0,4 eV, @o Gmpoverishment) =0,1...0,5 eV 


For lowering of temperature the efficiency of disperse particles as « snares» for electrons increase. For choosing 


nr 
at fig. 1 parameters the relative impoverishment of matrix 1— int) practically do not depend from € and make 
Ng 


up for temperatures 800, 1200, 1500 K conformably 76, 62, 54 %. 

b) Matrix enrichment (go < 0). 

The equation (4) for @o < 0 is essentially nonlinear, so solutions have been done numerical works. The 
calculative algorithm have been described in detail in [7]. For illustration at fig. 2 we bring calculations the 


distribution of electron density (the solution of equation (4)). 


Figure 2. Density of free carriers vs relative distance from boundary «matrix-particle» (Q = 1%, «= 10). 1-5 
nm, 2 - 10 nm, 3 - 20 nm, 4 - 50 nm. 

It’s shown that around metal particle "fur coat" from free carriers is forming which relative size depend 
essentially from ro. So if size of enrich zone R with density n = 500 no for particles with radius ro = 50 nm is R= 


1,5ro, then for ultradisperse particles with ro = 5 nm R = 4,5ro0. 


3. The Influence Of Particles' Ensemble On Electron Density In SOFC Materials. 


Lets go over to ensemble of disperse particles now. Disperse system, consisting from particles' ensemble in a 
matrix with extent enough areas of enrichment (impoverishment) and high particles' density such to produce area's 
interlinks, will be possesses heightened (low) density of free carriers in comparison with initial material and so 


heightened (low) matrix conduction. Properties of such system will be define in first approach by totality of single 


| 3 
particle solutions, if define radius R «influence area» of particle as R = 3 ZN , where N - particles’ density in 
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volume unit of matrix. In this case totality of influence areas will exhaust total volume. Selfco-ordinated account of 


«influence areas», taking into account particles sizes distribution, for equations, similar to (4), given in [8]. 
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Lets introduce Q volume fraction of dispersion particles, Q= 3 = 3 


small as possible to preserve basic matrix properties, for instance, ionic conductivity, mechanical properties etc. 


. Evidently, Q must be as 


We show upper limit Q = 30...40%, which, as well known, correspond to percolation threshold (under such 
conditions it is possible formation of through bridges from dispersion particles). This limit restricts applicability of 


examining model with isolated particles. Therefore volume fraction 1-10% was used in the calculations. 


Now we use obtained calculation data for prognosis of properties of dispersive compositions conformably to 
SOFC materials. The values of work of exit of electrons are shown in table | [9]. One can see from comparison of 
values y for some materials, their difference @o corresponds with calculation data. 

Table 1 


Values of works of exit of electrons y some materials[9] 


Work of exit, eV Throw of values, eV 


a 


Evidently, discussed in section 2a model can be utilised to the materials of solute electrolyte. As follows from 
table 1, pare parameters ZrO2-W are corresponding to calculation data (Fig.1). As it is follows from Fig.1 based on 
the calculation data Q = 1 %. Thus, in accordance with proposed model, introduction 1 % dispersive particles W of 
size ~10 nm in ZrO can decrease electron density approximately at one order. The second significant consequence 
of obtained data is effect increasing for low temperatures. This fact gives the possibility to change correlation 
between ionic and electronic conductivity to the side of ionic for low temperatures and therefore decrease ECG 
working temperature. 

Go over to discussion of case b) - matrix enrichment by free carriers. Using obtained data, one can estimate 
particle's density N, that required for zone's interlink and their volume fraction Q. When ro = 50 nm then Q=0,3 , 
when rp = 5 nm Q=0,01. 

Thus, when particle's size ro = 50 nm one can increase free carriers’ density in total volume on approximately 
one order. If ro = 5 nm that density increasing is possible long since Q = 0,01. When Q is a great value we can wait 
more substantial carrier's density increasing in nanocomposite material. Electron density change in rare-earth 
elements’ cobaltites, manganites, chromites [10] can be related with predicted enrichment of volume by free 
carriers. At last, using of obtained data allow examine as an envisaging further development electrode materials 


other class of substances. 


4. Irradiation-Plasma Methods Of Forming Nanostructures. 


Considering of the real ultradisperse structures conformably to SOFC materials one have to take into account 


some additional factors that become complicated problem, in particular, dispersion of the dimension particles, 


uniformity of their distribution matrix, thermal stability etc. Apparently the state of the boundary «matrix-particle» 
is imagine factor. The existing defect structure of the interface boundary is determine dispersion of emission data 
[9] and serves as base of different model of nanostructures [4]. It makes high demands to the technology of 
nanocomposite's manufactory. Among the most promising methods it is necessary to distinguish processes that are 
1) do not introduce (or introduced minimum amount) of impurities and 2) realised the disperse particles forming at 
high supersaturation of point defects. First of all one must consider vacuum and in situ methods. 

The structure properties of the two-phase materials Me - ZrO2 (Me = Ta,W,Re,Nb,Mo etc.) which received by 
vacuum CVD and PECVD methods have been investigated [11-13]. As electron microscopy showed the structure 
these materials represent two-phase system: matrix with uniformly distributed 5-10 nm. The implantation methods 
are very interesting to manufacturing of nanostructures in situ. The effect of forming of ultradisperse fraction under 
low temperature irradiation in multicomponent systems has been predicted [14]. The theoretical model has been 
experimentally checked [15]: the original particles have been solved under ionic bombardment of the nimonic PE- 


16 as new fine particles with the sizes 10 nm yielded. 


5. Conclusions 


5.1. Theoretical analysis was performed of electron density distribution n = n(r) around single particle with work 
of exit of electron differ from matrix. The dependencies n(r) are obtained at different values of parameters: 
temperature, @o etc. The concentration of electrons in matrix near particles are a) decrease for @o > 0 and b) 
increase for Qo < 0. 

5.2. The influence of particles’ ensemble on the changes in electron density in solid electrolytes has been studied. 
It is shown that addition of ~ 1% ultradisperse (~ 10 nm) particles with @o = 0,3 eV decrease electron concentration 
in matrix at 1500 K approximately to one order. The effect becomes intensive at lower temperatures. It is possible 
to suggest such nanocomposite materials for solid electrolytes with lower working temperature. 

5.3. The distribution of electron density n in nanocomposites with @o < 0 have been considered numerically. As 
calculations shown, addition 1 % fine particles with radius 5 nm increased n more those two orders. This effect 
gives opportunity to expanse the types of materials for electrodes (first of all, air electrode) SOFC. 


5.4. The irradiation-plasma methods for nanocomposites manufacturing have been examined. 
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